arXiv:astro-ph/0512487v2 23 Jun 2007 


Semilocal Cosmic String Networks 


Ana Achiicarro, 1,2 Petja Salmi, 1,3 and Jon Urrestilla 3,4 

1 Instituut-Lorentz, Universiteit Leiden , P. O. Box 9506, 2300 RA Leiden, The Netherlands 
2 Department of Theoretical Physics, University of the Basque Country UPV-EHU, 48080 Bilbao, Spain 
3 Department of Physics & Astronomy, University of Sussex, Falmer, Brighton BN1 9QH, UK 
*Institute of Cosmology, Department of Physics and Astronomy, Tufts University, Medford, MA 02155, USA 

(Dated: February 5, 2008) 


We report on a large scale numerical study of networks of semilocal cosmic strings in flat space 
in the parameter regime in which they are perturbatively stable. We find a population of segments 
with an exponential length distribution and indications of a scaling network without significant 
loop formation. Very deep in the stability regime strings of superhorizon size grow rapidly and 
“percolate” through the box. We believe these should lead at late times to a population of infinite 
strings similar to topologically stable strings. However, the strings are very light; scalar gradients 
dominate the energy density and the network has thus a global texture-like signature. As a result, 
the observational constraints, at least from the temperature power spectrum of the CMB, on models 
predicting semilocal strings, should be closer to those on global textures or monopoles, rather than 
on topologically stable gauged cosmic strings. 


I. INTRODUCTION 

The improvement in observational data during the 
last years, especially on Cosmic Microwave Background 
(CMB) radiation, has shown that topologically stable 
cosmic strings cannot be the dominant component creat¬ 
ing the primordial fluctuations that seed large scale struc¬ 
ture as we see it in the Universe today. Strings could be 
a subdominant component [l[ but there is as yet no ev¬ 
idence for their existence (see [2:] for a candidate string 
lensing event that turned out to be two galaxies). On 
the other hand, from the theoretical point of view, the 
formation of topological defects via the Kibble mecha¬ 
nism p] is inevitable in theories that allow stable defects. 
The formation of cosmic strings after inflation appears 
to be generic both in supersymmetric grand unified the¬ 
ories breaking down to the Standard Model ill and in 
brane inflation models within the framework of funda¬ 
mental superstrings @]. Therefore, models that do not 
predict topologically stable defects at the end of infla¬ 
tion are somewhat appealing. Recent work {?)] suggests 
that a particular type of non-topological defects, called 
semilocal strings, could form in the early universe both 
after U-term inflation f§] and after brane inflation with 
L>3/D7-branes @,E- 

The semilocal model El is a minimal extension of 
the Abelian Higgs model by a global 517(2) symmetry. 
With two equally charged Higgs fields and only one 17(1) 
gauge held, there are not enough gauge degrees of free¬ 
dom to cancel all scalar gradients, even outside the string 
cores, unless the scalar windings are correlated. Semilo¬ 
cal (SL) strings are stable for sufficiently large gauge cou¬ 
pling EL E 1 but due to their non-topological character 
they appear after the phase transition as open segments. 
These are closely related to electroweak dumbbell con¬ 
figurations E!, H, EB| but whose ends have long-range 
interactions that resemble those of global monopoles [1% 
with a force roughly independent of distance. Therefore 
the network of SL strings has several features that are 


not present in networks of topologically stable Abrikosov- 
Nielsen-Olesen (ANO) strings [17} and their cosmologi¬ 
cal evolution could be quite different EE- Semilocal 
string segments can either contract and eventually disap¬ 
pear or grow to join a nearby segment and form a longer 
string. Closed loops can be formed by intercommuta¬ 
tions, and also if the two ends of a segment join. The 
obvious question is if at some point this joining of seg¬ 
ments can form infinitely long strings. If so, the networks 
present characteristics intermediate between topological 
ANO strings and global defects - both of which show 
scaling behaviour in an expanding universe E - and we 
would like to understand their cosmological impact. 


The purpose here is to report on a large scale numerical 
study of SL string networks. The non-topological nature 
of the SL strings does not permit a Nambu-Goto ap¬ 
proach and our results are based on classical field theory 
simulations. Our basic results are twofold. We find that 
SL strings will, in the right parameter range, form a net¬ 
work that shares features with a network of ANO strings 
due to the existence of fast growing, extremely long per¬ 
colating strings, which we believe would lead to an in¬ 
finite string component at late times (we use the word 
“percolating” in a non-technical sense, to mean that the 
strings span the simulation box). On the other hand, the 
energy density is dominated by scalar gradients, like in 
models for global defects, even in the case of dense string 
networks deep in the stable regime. As the energy den¬ 
sity plays a crucial role for the observable cosmological 
consequences, like the effect on the temperature power 
spectrum of the CMB, the signature of a SL string net¬ 
work is thus likely to be very similar to that of global 
defects. 
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II. SEMILOCAL MODEL 


The action in dimensionless units is given by 


S = <F x 


D,M 2 - -F U „F^- 


0 (^- 1 )' 


,( 1 ) 


where D M — iA^ and d> T = {(f> 1 ,(^ 2 ), is a complex 

SU{ 2) doublet. The only parameter of the model is the 
ratio of the scalar and gauge masses, (3 = TO 2 calar /TO gauge 
(more generally, (3 oc X/q 2 , where q is the gauge coupling 
and A the quartic coupling). 

ANO flux tubes embbeded in the SL model (e.g. with 
<f >2 = 0 ) are classical solutions of (P) but they are not 
topologically stable because the vacuum manifold (the 
three-sphere |0i| 2 + \(j>2 \ 2 = 1) is simply connected. A 
string with <f> 2 = 0 can decay into configurations with 
<f>2 7 ^ 0 and stability is a dynamical question. It was 
shown in [llj, [l 2 | that infinite straight vortices are sta¬ 
ble if /3 < 1 (neutrally stable at (3 = 1) for any winding 
number. For (3 > 1 the bare embedded ANO strings 
with unit winding are unstable toward spreading of the 
magnetic flux. (A new class of solutions carrying per¬ 
sistent currents was found [ 20 | whose stability is under 
investigation). 

The string core carries U( 1) magnetic flux. We use 
this magnetic field to identify the strings in numerical 
simulations, as observing the zeros of scalar fields does 
not provide an unambiguous detection method for non- 
topological defects [l9} . A SL string can be continuously 
deformed into “skyrmion-like” configurations; an excited 
string may have no scalar zeroes at the core but neverthe¬ 
less it is a linear concentration of energy and for cosmo¬ 
logical applications it should be considered a string-like 
defect. Here we follow the strategy of [H, (23|: we sep¬ 
arately obtain numerically the profile of the scalar and 
magnetic field of the straight ANO string for several val¬ 
ues of (3 and use this information to analyze the data 
from simulations. If the magnetic field strength at a 
given point exceeds a certain fraction / of the maximum 
magnetic field of the corresponding ANO string core, the 
point is considered to be part of a SL string. The string 
core fattens with decreasing values of (3 , therefore we de¬ 
fine the length of a string to be the number of lattice 
points belonging to the string divided by the cross sec¬ 
tion of the corresponding ANO string. We will discard 
blobs whose length is less than a cutoff (set to be order 
of the width of a string). 


III. SIMULATIONS 

Previous numerical studies of SL strings have tested 
their stability [ 2 l[ and formation rate in two [ 22 } and 
three [23j dimensions, showing that formation rates can 
reach a third of the rate of topologically stable strings at 
the lowest values of (3 simulated (which was (3 = 0.05). 
Based on these studies, we work under the assumption 


that the SL string network forms and we utilise present 
high performance computing technology to achieve much 
longer dynamical range to study the evolution of the net¬ 
work. 

The discretized field equations were derived from (CD 
using techniques from Hamiltonian lattice gauge theories 
[24l | and solved using the staggered leapfrog algorithm. 
Damping terms r]cj) and r/A^ are introduced into the scalar 
and gauge held equations, respectively. We set the lattice 
spacing to be Ax =1.0 (the string cores are at least three 
points wide) and the time step At = 0.2. 

During the simulation total energy densities are cal¬ 
culated and strings are monitored via their magnetic en¬ 
ergy. We study two main string observables as a function 
of time: the distribution of strings as a function of length 
n(l, t) and a length scale £(f) = a, JV/L(t) where V is the 
volume and L(t ) = / 0 °° n(l, t)ldl is the total string length. 
Linear growth of £ as a function of time is referred to scal¬ 
ing of the string network. With periodic boundary condi¬ 
tions and our choice of lattice spacing the time required 
for two signals emitted from the same point and travel¬ 
ling in opposite directions to interfere with each other is 
half of the side of the grid. This sets the maximum time 
the simulations can probe reliably. If a string spanning 
the box forms before that time, from the point of view 
of an observer in the box it would appear coming from 
superhorizon scales. 

We are interested in the regime (3 < 1. Fattening of 
the scalar core sets a lower bound \f]3 = 0.2 that we 
can examine numerically; for smaller values of [3 scalar 
cores of different strings start to overlap in the initial 
configurations used. We tested several choices of /, the 
threshold of the magnetic field that decides if a point 
belongs to a string. For 0.2 < / < 0.4 we do not observe 
any significant discrepancy in the length scale £ of the 
network. 

Based on earlier works [U, m we do not expect the 
formation of the network to be very sensitive upon the 
initial conditions. Two different strategies for setting the 
initial conditions were considered: 

a) One strategy consists of setting all fields (scalars 
and gauge) and gauge momenta to zero but giving the 
scalar field momenta some velocity obtained from a set 
of Gaussian random variables with zero mean. The con¬ 
figuration of uncorrelated velocities is then smoothed by 
averaging the velocities over nearest neighbours s times. 

b) The second strategy is closer to the Vachaspati- 
Vilenkin simulations [25l | where all held momenta an 
gauge fields are set to zero, but random phases are 
asigned to the scalar held, which in turn are smoothed 
out by averaging over neighbouring points s times. 

Note also that both types of initial configurations sat¬ 
isfy Gauss’ law. A set of simulations was performed in 
simulation boxes of size 384 3 using both kinds of initial 
conditions varying the parameter s from 10 to 40. We 
monitored the length scale £ and found that the simula¬ 
tions did not show qualitatively different behaviour, and 
quantitatively, the deviations lay within the statistical 
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FIG. 1: A snapshot of the simulation box (512 3 ) for y 7 /? = 0.2, 
rj = 0.05 at t = 250. 


errors obtained by repeating the simulations with several 
(of order ten) different initial configurations. Moreover, 
case b) has typically enough energy so that the fields 
climb up the potential and then the system goes to a 
state resembling that of case a). 

After some transient time, the system forms a network 
of strings, and the network reaches a scaling regime. This 
happens for all the initial conditions studied above, the 
only difference being the time it takes for the simula¬ 
tion to go through the transient regime. Therefore, we 
chose the most favourable way of creating initial condi¬ 
tions from the point of view of reaching quickly a scaling 
regime. As the dynamical range of simulations is in any 
case fairly limited, we seek to minimize the initial tran¬ 
sient time. We do not expect to capture all the physical 
phenomena in the phase transition with this treatment, 
but neither is the main focus of this study the precise 
description of the dynamics at the phase transition, but 
rather the evolution and properties of the network, if 
formed and persistent, after the phase transition. For 
instance, the initial conditions used are not thermal, but 
this is the case for e.g. tachyonic preheating after infla¬ 
tion [26j . 

In the early universe there would be several sources 
of dissipation for the string networks, most notably the 
Hubble damping that decreases inversely proportional to 
time. In the case of our simulations there is the need 
to set the damping in such a way that the subsequent 
evolution of the network proceeds fast but without the 
fields overshooting. Several choices for r) were tested and 
for the data shown it is either a constant (rj = 0.05), or 
inversely proportional to time ( r\ = 6/f). During the first 
time steps (t < 30) the damping is set to rj = 0.2 to go 
through the initial transient period fast and lose energy 
efficiently. 


IV. RESULTS 

Simulations in boxes size 512 3 were performed for a set 
of ten (y/d > 0.3) or fifteen (y//3 < 0.3) initial configu¬ 
rations for each value of y 7 /? ranging from 0.6 to 0.2 (we 
also verified the disappearance of the string network at 
the Bogomol’nyi limit, (3 = 1). We used initial random 
velocities with smoothing 20 times (s = 20), the frac¬ 
tion / = 0.25 for string identification and damping term 
either constant (rj = 0.05) or time dependent (77 = 6/f). 

After the initial transient (f ~ 30) a network of string 
segments is formed. The subsequent evolution of the seg¬ 
ments depends strongly on the value of /?, but as a gen¬ 
eral rule strings are longer and their density increases as 
/3 decreases. Intercommutations of strings are rare. 

For y/d = 0.6 the maximum string length is around 
half the size of the simulation box. In principle, two 
strings of that length joining together could form a string 
spanning the box size, but we have not seen any such 
event, because the network is very sparse. For y/3 = 0.4 
we have the first event of a string that spans the simula¬ 
tion box at time = 250. Decreasing (3 even more leads to 
longer strings and a more dense network, for y/d = 0.3 
approximately one third of the initial configurations lead 
to the formation of at least one percolating string; for 
v/d = o.2 more than a half. Figure [T] shows a snapshot 
of the simulation box for y/d = 0.2, where there are sev¬ 
eral strings spanning the size of the simulation box. For 
y/d <0.25 we find some extremely long strings. Analy¬ 
sis of the length distribution shows a qualitative change 
around this value of y/d ss 0.25. 

For y/d > 0.25 the segments lie in an exponential dis¬ 
tribution. Nevertheless, for y/d < 0.25 we find that be¬ 
sides the population of segments in the exponential dis¬ 
tribution, we find extremely long strings (Z/£ = 19 in 
Figure [2j which are clearly not part of the exponential. 
These long strings percolate through the box and are the 
precursors of infinite strings at late times. 

The value of (3 influences the network in two different 
ways. First, the initial string density (after discarding 
the transient) grows with decreasing [3, see figure [3] Sec¬ 
ond, the accretion of disperse magnetic flux into tubes is 
much more efficient for low (3. To try to disentangle the 
two effects we have tested configurations with different 
initial defect densities for low (3 and found that they all 
quickly approach the same £(f), suggesting that flux ac¬ 
cretion is the dominant effect that determines the scaling 
behaviour. 

When the string network becomes more dense, there 
are also occasionally loops forming, but approximately 
only every third simulation shows a loop for y/d < 0.4. 
The reason might be simply an observational effect, loops 
disappear fast and are not likely to be seen in snapshots. 
Moreover, we have not seen any intercommutation in our 
simulations. On the other hand field theory simulations 
do not tend to show many loops for ANO string networks 
either [27|. We do not observe a population of small 
loops and such a tiny number of loops suggests that they 
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P=0.04 n=0.05 



FIG. 2: The string length distribution averaged over 15 sim¬ 
ulations for y/j3 = 0.2, 7] = 0.05, at time = 150 (shaded) and 
250 (thick line). Up to l/£ ~ 10 — 13 the segments lie on an 
exponential distribution. The other points are extremely long 
strings not belonging to the exponential. Those cases are only 
found for \fp < 0.25. 

100 
90 
80 
70 
60 

sXP 

50 
40 
30 
20 
10 

50 100 150 200 250 

Time 

FIG. 3: £(t) in a simulation box in a single run with y/]3 = 
0.6, 0.5, 0.4, 0.3, 0.2 (from top to bottom), r/ ~ i -1 . 


cannot have a major effect in the network. 

Figure [3] shows the observable £ as a function of time 
in simulations where the damping term was set to be 
r/ = 6/1. The slow change in time for small /3 proves that 
the string network can persist for a considerably long 
time and goodness of linear fit to data indicates a scaling 
network within the range we can probe. 

Finally, we monitored the total energy densities in the 
simulation box. In all cases studied here the scalar gradi¬ 
ents dominate the energy density. This is not surprising 
in the Bogomol’nyi limit, since the string network disap¬ 
pears, but it is the case also even at lowest values of (3 
we were able to probe. The approximate ratios of energy 
densities ( kinetic : potential : magnetic : scalar gradi¬ 
ents ) are (2:1:1: 15) for yj~j3 = 0.5 and (1 : 1 : 1 : 7) 
for yfP = 0.2. 


V. DISCUSSION 

Consider the exponential distribution of short seg¬ 
ments. At low [3 the string segments are very light com¬ 
pared to the (global) monopoles at the ends, which are 
strongly interacting and playing substantial role in the 
evolution of the segments. This is just the opposite of 
the usual situation encountered in cosmological scenar¬ 
ios in which the magnetic monopole problem is solved 
by a subsequent phase transition in which the magnetic 
flux of the (local) monopoles is confined to strings. In 
this latter case the monopoles are very light and hardly 
contribute to the dynamics, which are dominated by the 
string tension. 

If a semilocal segment is much shorter than the average 
interstring separation it is very likely that the segment 
will shrink and the monopole and antimonopole at the 
ends will annihilate each other. In this case one naively 
expects similar length distributions whether the strings 
are heavy or light. However, for long segments the dif¬ 
ference becomes apparent. In particular, the breaking 
of a semilocal string into a monopole-antimonopole pair 
is strongly suppressed [28(. Tunneling is forbidden be¬ 
cause there is no configuration that the straight string 
can tunnel to, as the string configuration has lower energy 
than the broken configuration. In this respect, semilocal 
strings at low (3 resemble topological strings. The break¬ 
up instability is the reason why in the usual cases of 
(light) monopoles connected by strings, the length dis¬ 
tribution stays exponentially supressed throughout. By 
contrast, in the semilocal case there is not such a sup¬ 
pression of long strings. Long SL strings, once formed, 
are likely to continue growing by joining. 

Thus it would be tempting to argue that, deep in the 
stability regime, a SL string network could be considered 
effectively as a collection of global monopoles joined by 
comparatively light strings. If this picture is valid, the 
scaling in cr-model studies of global monopoles in an ex¬ 
panding universe pf)l . [30| provides further support for the 
scaling behaviour we observed within the limited range 
we were able to probe. As time progresses, for sufficiently 
low (3 we expect the appearance of a population of infi¬ 
nite strings with typical ANO behaviour, loop produc¬ 
tion, etc. but at a lower density than in the topological 
case starting from the same initial state. 

VI. CONCLUSIONS 

In this paper we studied the evolution of semilocal 
string networks deep in the stability regime using high 
performance computing to carry out numerical field the¬ 
ory simulations in flat space. Our main results are: 

• A population of segments with an initially exponen¬ 
tial length distribution. Short segments vanish rapidly 
and the mean string length grows. 

• Appearance of superhorizon strings (much larger 
than the box size) percolating through the box deep in 
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the stability regime, y(9 < 0.3. For yffi < 0.25 these fast 
growing strings fall outside the exponential distribution. 

• No significant population of small loops. In dense 
networks \JJ3 < 0.4 about one third of the simulations 
show a loop of a larger size. 

• The string network can persist for a long time and 
for low 0 the results are consistent with linear scaling. 

• The energy density is dominated by scalar gradients. 

These findings suggest that SL networks may have a 

cosmological signature, in particular in the CMB and 
large scale structure, closer to that of global defects 
(monopoles or textures) than to topological strings. Cur¬ 
rent bounds [53] on global textures in CMB allow about 
13% of the signal to come from defects at £ = 10. 

The relevance of these results to the cosmological sce¬ 
narios in 00 would require further study in expanding 
backgrounds (see also [l8|) and the development of an¬ 
alytic characterizations of the network to verify the ap- 
perance of an infinite string population at late times, as 
was done in ref. [32} for the case of dense networks of 
cosmic string loops. Apart from the expansion, a com¬ 
plete quantum field theoretical study of the conditions at 
formation is still lacking, as pointed out in Q. 

If it is confirmed that the observational signal of 
the semilocal network is similar to global textures or 
monopoles in the temperature power spectrum [33], there 


exists a potential problem for distinguishing between 
models via future data. A dedicated study is needed to 
find out if there is any difference e.g. in polarisation of 
vector modes, or in the possibility of gravitational waves 
radiating from oscillating loops [17} ■ 
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